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O Abstract 

o 

^ We present here an approach to a computation of C(2) by changing 

i-H variables in the double integral using hyperbolic trig functions. We also 

|JZ ^-Pply this approach to present Cl'^); when n > 2, as a definite improper 

Qs integral of single variable. 
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The Riemann zeta function is defined as the series 



' ' 1 Introduction 

> 

cn 

en C{n) 

O for any integer n > 2. Three centuries ago Euler found that ({2) = vr^/G, which is 

an irrational number. Exact value of ({3) is still unknown though it was proved 
by Apery in 1979 that ^(3) was also irrational (see f5]). Values of Cl^)? when n 
is even, are known and can be written in terms of Bernoulli numbers. We refer 
the interested reader to chapter 19 of p] for a "perfect" proof of the formula 
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n=l 



^2'= {2k)\ 



Notice that ({n) can be written as the following multi- variable integral 



1 rl 



({n) = I ■I dxidx2 . ■ ■ dxr, 

'o -'- ~ 2:1X2 ■ ■ ■ Xn 



1 



Indeed, each integral is improper at both ends and since the geometric series 
^ x'^ converges uniformly on the interval |a;| < i?, Vi? G (0, 1) we can write 

9>0 



^^{xiX2---Xn)'' 



1 - a;ia;2 ■ ■ ■ Xn ^^ 

then interchange summation with integration, and then integrate {xiX2- ■ ■ XnY 
for each q. Using the identities 

11 2 , 1 1 2x1/ 

and 



1 — xy 1 + xy 1 — x'^y'^ 1 — xy 1 + xy 1 — x'^y'^ 

and a simple change of variables one can easily see that 

/ / -. dxdy = - / ^^ dxdy 

Jo Jo ^-xy 3 Jo Jo 1 - X y 

by further generalizing this idea one comes along the following, 

^n /•! /■! ^ 



CW = ^^ T / • • • / n dXi...dXn. 

^ Wo -^0 l-ria;? 

1=1 



Notice that (1, 1) is the only point in the square [0, 1] x [0, 1], which makes the 
integrand 1/(1 — x^y^) singular. If we take another point on the graph of 1 = x^y^, 
say (a, 1/a) with a G (0, oo), then it follows easily (see lemma 1 below) that 

l/a r-a i C^ f^ A 

Jo 



„ „ dxdy = TT^ dxdy. 

1 - x^r Jo Jo 1 - X y 



This result motivates the following definition 

Definition 1. For any point (oi, 02, ... , dn-i) ^ M"^"*^ such that ai G (0, +00), Vi G 
{1, . . . , n — 1} we define 

/„(ai,...,a„_i) = / ... / / dxidx2. . .dxn- 

Jo Jo Jo i_Yl xl 

Lemma 1. For any ai G (0, +00), we have /n(ai, ..., fln-i) = In{^, 1? • • • > !)■ 



Proof. Simply observe that by using the change of variables Xi = aiUi for all i G 
[1, . . . , n], where an = l/(aia2 ■ ■ ■ On-i), the Jacobian equals 1, and the integrand 
is unchanged. D 

In this article we investigate ({n) following Beukers, Calabi and Kolk (see [2]), 
who used the change of variables 

sinfti) , sin(f) , f^ f^ 1 , , 

X = -— and y = -— to evaluate / / — — ■ axay. 

cos(t;) cos(m) Jo Jo 1 - ^ V 

Such proof of the identity ({2) = tt'^/6 may also be found in chapter 6 of [1] and 
in papers of Elkies [3] and Kalman ^4j. Let us also mention here that Kalman's 
paper, in addition to a few other proofs of the identity, contains some history of 
the problem together with an extensive reference list. 

Here we will be changing variables too, but in the integrals /n(ai, •••, o-n-i) and 
using the hyperbolic trig functions sinh and cosh instead of sin and cos. Such 
a change of variables was considered independently of us by Silagadze and the 
reader will find his results in [B]- 

Acknowledgement: The authors would like to thank the referee who drew 
our attention to Kalman's paper [3] and has made a few useful suggestions that 
improved the exposition. 



2 Hyperbolic Change of Variables 



First observe that the change of variables 

sinf-Uj) , ,. ^^ , , , 

Xi = T^^ Vz G N mod (n) 

COS{Ui+i) 

reduces the integrand in /„(!,..., 1) to 1 only when n is even. The region of 
integration $„ = [{xi, . . . , x„) G M" : < Xi, . . . , x^ < 1] becomes the one-to-one 
image of the n-dimensional polytop (note Un+i = Mi) 

71" 

n„ := [(mi,M2,...,m„) G M" : Mj > 0,Ui + Ui+i < -,1 < i < n]. 

We suggest here a different change of variables that will produce an integrand 
of 1 for all values of n in /„(ai, ...,0^-1). But first we define the corresponding 
region. 



Definition 2. For any point (oi, 02, ... , On-i) G M" ^ such that ai G (0, +00), Vi G 
{1, . . . , n — 1} we define 

<l>„(ai,a2,...,a„_i) := [(xi, . . . ,x„) G M" | < Xj < a^, Vi G {l,...,n}], 
where an = l/(ai ■ 02 ■ . . . ■ a„_i). 
Lemma 2. The change in variables 

sinh(ui) 



Xi 



cosh(u 



Vi G N mod (n) 



i+ly 



reduces the integrand of In{ai, ...,a„_i) to 1 for all values of n > 2. It also gives a 
one-to-one diff'erentiable map between the region $„(ai, 02, ... , a„_i) and the set 
Tn C M" described by the following n inequalities: 

< Ui < arcsinh(ai ■ cosh(ui+i)), Vi G N mod (n). 



Proof. The inequalities for r„ follow trivially from the corresponding inequalities 
< Xi < Qi and the facts that cosh(x) > and arcsinh(x) is increasing every- 
where. Injectivity and smoothness of the map may be proven by writing down 
formulas, which express each Ui in terms of all Xj. For example, here are the 
corresponding formulas for the set T^: 



Ui = arcsmh Xi ■ ^ / ^ r, ^ ■ , 

1 - xf x^x^ 



i G N (mod 3). 



The Jacobian is the determinant of the matrix 

/ cosh(Mi) — sinh(Mi) sinh(M2) 



A 



cosh(M2) 





- sinh(u„) sinh(«i) 
cosh (mi) 



cosh (M2) 
cosh(«2) 
cosh(u3) 







' sinh(M2) sinh(M3) 
cosh (M3) 



\ 





cosh(M„) . 
cosh(Mi) / 



To compute this determinant we observe that the first column expansion re- 
duces the computation to two determinants of the upper and lower triangular 
matrices. This results in the formula, where the first term comes from the upper 
triangular matrix and the second from the lower triangular matrix (recall that 

Un+l = ui) : 



Det(A) = fr ''f'}''^\^ + (_i)n-i.-[T -^inh(. )sinh(..,,) ^ ^_-rT^^^^^ 

t\cOsh{Ui+i) f\ COsh'(M,+i) -LJ- 



'%). 



i=l 



When using the above change in variables the denominator of the integrand 



n n 

,2/ 



1 — J]^ Xj becomes 1 — H tanh (mj), which we just proved to be the Jacobian. D 



3 Computations of ({2) 



We begin with C(2), which is a rational multiple of /2(1)- Lemma 1 implies that 
it's enough to compute 



/2(a) = / / :^-^dxdy for arbitrary a > 0. 

'0 Jo 



1 — x'^y'^ 



We now preform the following change in variables 

sinh(u) sinh(v) 



^ = — T^T' y 



cosh(v) ' cosh(u) 

As we proved above, our integrand reduces to 1 and all we must do is worry 
about the limits. If a; = then clearly u = 0, the same is true for y and v. 
If X = a then a ■ cosh(v) = sinh(u) so f = arcosh ^^'" ^^' ) and ii y = - then 

(1/a) ■ cosh(u) = sinh(v) so f = arcsinh( '^°'^ ^"^ ) thus describing our region of 
integration (see Figure 1). We then write the integral /2(a) as follows 

arcsinh(a) oo 

arcsinh ( | du+ / arcsinh ( 1 — arcosh ( I du. 



arcsinh(a) 
r<arcsinh(a) . , / cosh(u) ^ 



Lemma 3. lim P'''^^™ arcsinh ( '^"'^ ^"M du = 



a->0 



Proof. If we let cosh(arcsinh(z)) = Q then Q = \/l + z'^. Therefore 



/cosh(arcsinh(a))\ , / / 1 

arcsmh | I = arcsmh \/ ^ + 1 




arcsinh(a) 



Figure 1: The set Ta C R^, Va > 0. 



Since arcsinh(cosh(u)/a) is concave up, we can take area of the rectangle with 
vertices at (0,0), (arcsinh(a), 0), and (arcsinh(a), arcsinh(cosh(arcsinh(a))/a) as 
an overestimate of the integral, that is 



, , , 1 \ /•— i'lMa) /cosh(u) , , 

arcsinh(a) ■ arcsmh I a/ ^7 + 1 I ^ / arcsmh | | du > 



a I .10 

Then by applying L'hospital's rule one can deduce 



lini(arcsinh(a) ■ arcsinh 



a-s>0 



1 



in = o. 



n 



Now, since /2(a) = hi^), Va > 0, we conclude that hi^) = liniJ2(a), and 

a— !>0 

therefore we have 



/2(1) = hm 

a— >0 



arcsinh ( 1 — arcosh ( 1 du. 



arcsinh(a) 



Since 



and 



we get 



, /cosh(x)\ , / cosh(x) /cosh^fx) 
arcsmh ( ^^ I = In ( ^^ + d ^ + 1 



, ,'sinh(x)\ , / sinh(x) /sinh^fx) 
arcosh ^—- = In ^—- + \ ^^^ - 1 



a2 



^ / cosh(x) / cosh2(x) j^ 

Jsfl) = lim / In -^ ^ "' I dx 

arcsinh(a) y smh(x) _|_ / sinh2(x) _ ^ 



a-5>0 

a 



which, after taking the hmit as a — )■ gives 



«^)-rKS|y'^- 



Using integration by parts -u = In ( ^°^^K^ j and v = dx one obtains the formula 

'cosh(x) 



Io(l) = x\n 



sinh( 



X 



sinh(2x 



2x 

dx. 



By examining the hmits of the first half of the formula as x goes to and oo we 
are left with only the integral 

f°° 2x 
Jo smh(2x) 

By applying the change in variables u = 2x our formula becomes 

1 f°° u 
^2(1) = 7: . , . . du. 

2 Jq smh(u) 

Now we use the method of differentiation under the integral sign and consider 

the function 

^, , 1 f°° arctanhfatanhfx)) , 
Fa = - / -— dx. 



2 7o sinh(xj 

One should consider the function F at the points a = 1 and a = 0. F{1) is 
clearly the integral we are trying to find and -F(O) is 0. Thus by differentiating 
under the integral with respect to alpha, plus some algebra we obtain 

w/ X ./ X 1 P cosh(x) 

F (a) = f{a) = - -. ^7 —T——dx. 

^ ' ■'^ ' 2 Jo 1 + (1 - a2)sinh2(x) 

7 



Then by preforming the change in variables u = y/l — a^ ■ sinh(x) the integral 

becomes 

1 /"°^ 1 

2V 1 - a^ Jo 1 + M^ 



which is simply 



arctanh(u) 



2y/T^ 



a^ 



n 



4y/T^ 



Ct^ 



Since we took the derivative with respect to a we must take the integral with 
respect to alpha so we have 

• 1 
f{a) da = F(l) - F(0) = F(l) - = F(l) 



^0 

which, as stated above is our goal. So 








HI) = f f\a) da = lf -=i= da = 
Jo ^ Jo V 1 - d^ 


— arcsinfa) 
4 


1 






and hence C(2) = | ■ ^ = vrV6. 









4 A formula for C(^)) n>2 



One could try to use similar approach to compute C(^)) n > 2, however the 
computations become a bit long. Instead, we present an elementary proof of 
the following theorem, which generalizes our formula for C{2) from the previous 
section. 



Theorem. Let n > 2 be a natural number. Then 

dxi...dxn = 7 7T7 ■ / ln"~^(coth(x)) dx. 



1 ^1 1 -1 ^oo 



JO i_j|3.2 (n-1)! Jo 

i=l 



Let us start with the following lemma, which can be easily proved by using 
induction on k, integration by parts and I'Hospital's rule. 

Lemma 4. 

Jo (2g+l)^'+i' 



Proof of the theorem. Applying the substitution z = tanh(x) to the integral 

^—^l ln-'(coth(x)) d. 
gives 

Since the integral is improper at both ends and the geometric series ^ z'^'^ con- 
verges uniformly on the interval \z\ < R, Vi? G (0, 1), the last integral equals 

— -^V(-l)""i- [ \n''~\z)z^''dz = by lemma 4 = V- ^^. 

Using the geometric series expansion one can easily show that we also have 
Jo -^0 1 - n a;2 ,>o ^2? + ^^ 



i=l 



D 



Corollary. For any integer n > 2, 

on ^co 



((n) = ^- • / In'^-Vcoth(x)) dx. 
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